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Abstract 

The positron anomaly recently reported by the cosmic-ray measurements 
can be explained by the decaying dark matter scenario, where it decays mainly 
into leptons with the lifetime of O(10 26 ) second. When the dark matter is a 
fermionic particle, the lifetime of this order is known to be obtained by a 
dimension 6 operator suppressed by the unification scale (~ 10 16 GeV), while 
such decay operators do not necessarily involve only leptons. In addition, the 
scenario would be spoiled if there exist lower-dimensional operators inducing 
the dark matter decay. We show in this letter that a single non-Abelian discrete 
symmetry such as A4 is possible to prohibit all such harmful (non-leptonically 
coupled and lower-dimensional) operators. Moreover, the dark matter decays 
into charged leptons in a flavor-blind fashion due to the non-Abelian flavor 
symmetry, which results in perfect agreements not only with the PAMELA 
data but also with the latest Fermi-LAT data reported very recently. We also 
discuss some relevance between the discrete symmetry and neutrino physics. 



1 Introduction 



The existence of non-baryonic dark matter, which accounts for about 23% of the 
energy density in the present universe, has been established thanks to the recent 
cosmological observations such as the WMAP (Wilkinson Microwave Anisotropy 
Probe) experiment pp. The detailed nature of the dark matter is, however, still un- 
revealed and is a great mystery not only in astrophysics and cosmology but also in 
particle physics. In order to detect and study the dark matter, various theoretical and 
experimental efforts have been devoted, and possible signals for the dark matter have 
recently been reported from the indirect detection measurements at the PAMELA (a 
Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics) [2J and 
Fermi-LAT (The Fermi Large Area Telescope) [31 H] experiments, where anomalous 
excesses of cosmic-ray positrons (electrons) have been found. Though it is under 
debate whether these anomalies are interpreted as dark matter signals, they have 
motivated many theoretical study to explore the nature of the dark matter. 

There are several types of scenarios to explain the cosmic-ray anomalies. Among 
these, we focus on the decaying dark matter scenario [51 El [7] where the dark matter 
is assumed to be unstable with the lifetime much longer than the age of the universe 
and its decay in the halo of our galaxy explains the anomalies. The observational 
data of positron (electron) excesses as well as the non-observation of anti-proton 
excesses in the cosmic ray [5] suggest that the dark matter mass should be on the 
TeV scale and it decays mainly into leptons with the lifetime of O(10 26 ) sec. An 
important question for this scenario is why the lifetime is so long, in other words, 
what is the origin of meta-stability of the dark matter. An attractive answer is that 
the meta-stability is derived from very high-energy physics such as Grand Unified 
Theory (GUT). When the dark matter is a TeV-scale fermionic particle, it could 
decay thorough a four- Fermi operator suppressed by the GUT scale A ~ 10 16 GeV 
and the width is estimated as F ~ (TeV) 5 /A 4 ~ 10~ 26 /sec. With this interesting 
relation between A and T, various explicit studies have been performed so far [5J. 

In the context of decaying dark matter, it seems however difficult to realize the 
main decay mode contains only leptons, not hadrons. In addition, there generally 
exists lower- dimensional operators inducing the dark matter decay, and then the 
estimation of T may be disturbed. A reasonable solution to these problems is to 
implement appropriate symmetry which forbids the rapid and/or non-leptonic decay 
of the dark matter. In this letter, we point out that the leptonically-decaying dark 
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matter is guaranteed with use of non-Abelian discrete symmetry acting on the gen- 
eration space. We focus on the A 4 flavor symmetry [H] as the simplest example. The 
dark matter is assumed to be a Majorana fermion that is singlet under A4 and the 
standard-model gauge groups. Identifying the effective decay operators and evaluat- 
ing the positron (electron) flux from the decay, we show that the A 4 invariance leads 
to a novel flavor pattern of the dark matter decay which well describes the cosmic-ray 
anomaly reported by the PAMELA collaboration. It also turns out that the total 
electron and positron flux is in perfect agreement with the latest Fermi-LAT data 
reported very recently [3]. We also discuss some relevance of discrete symmetry on 
the neutrino physics, i.e. the masses and generation mixing of neutrinos. 

2 Decaying dark matter and discrete symmetry 

In addition to the standard-model fields, a gauge-singlet fermion X is introduced as 
the dark matter (DM) particle. We assume that the baryon number is preserved at 
least at perturbative level. It then turns out that there exist various gauge-invariant 
operators up to dimension 6 [TU] which induce the DM decay: 



Dimensions 


DM decay operators 


4 


LH C X 


5 




6 


Lelx, H^HLH c X, (h c yd^h c e^x, 




QDLX, UQLX, LDQX, U^DE-fX, 




D^H C D^LX, D^D^H C LX, 




B^ u Lcr^H c X, W* v La^T a H c X 



Table 1: The decay operators of the gauge-singlet fermionic dark matter X up to dimension 
6. Here, L, E, Q, U, D, and H denote left-handed leptons, right-handed charged leptons, 
left-handed quarks, right-handed up-type quarks, right-handed down-type quarks, and 
higgs field, respectively (H c = eH*). On the other hand, B^, W®, and are the field 
strength tensor of hypercharge gauge boson, that of weak gauge boson, and the electroweak 
covariant derivative. 

This general operator analysis shows that the dark matter X can decay into not only 
leptons but also quarks, higgs, and gauge bosons at similar rates. Furthermore, a 
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quick decay of DM is induced if the dimension 4 Yukawa operator LH C X is allowed. 
One may try to impose an Abelian (continuous or discrete) symmetry to prohibit 
unwanted decay operators, but it does not work. The reason is the following: the 
abelian charges of L, E, and H are assigned to be q L , q E , and q H , respectively. 
The operator LHE should be invariant under the symmetry in order to have the 
masses of charged leptons, and the relation q L = q E + q H is hold. The invariance 
of LELX is also needed because this is the unique operator in Table [I] for the 
leptonic decay of dark matter, and leads to 2q L = q E + q x . These charge relations 
turn out to imply that q L + q H = q x and the operator LH C X necessarily becomes 
symmetry invariant. The discussion is unchanged even when the charged-lepton 
Yukawa coupling is generated from higher-dimensional effective operators, in which 
case, an unfavorable decay via LH C X is found to be suppressed by at most (electron 
mass) / (electroweak scale) and still leads to a short lifetime. Further, in Tabled], 
there are other leptonic decay operators such as H^HLH C X which do not contain 
hadrons. However they have the same property as LH C X with respect to the Abelian 
charge. 

In the following, we show that the desirable DM decay is guaranteed with use 
of non-Abelian discrete symmetry. Namely, non-Abelian symmetry allows us to 
prohibit the dangerous dimension 4 operator as well as other operators leading to 
non-leptonic DM decay, while keeping the operator LELX invariant. In this letter, 
we present a model with the discrete symmetry A 4 , though it is possible to construct 
different models with similar DM decay using other discrete symmetry. The A 4 group 
has one real triplet 3 and three independent singlet representations 1, 1', 1" [11] . 
The multiplication rules of these representations are as follows; 

3 © 3 = 3©3©1©1'© 1", 3 © 1' = 3 © 1" = 3, 

i'©i' = r", i" <g> i" = l', i'<g>i" = i. (2.1) 

One notice is that the multiplication of two 3's contains both 3 and real singlet 
1, and hence any products of more than two 3's can be invariant under the A4 
transformation. 

With this property of the A 4 symmetry, we consider the A 4 charge assignment 
given in Table [2j Remarkably, all the decay operators in Table [1] except LELX 
are forbidden due to this single symmetry, and the dark matter mainly decays into 
leptons. With the notation Lj = ((v e , e L ), (u^, fi L ), (v T , r L )) and E { = (e R ,ii R ,r R ), 
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Q U 


D 


L 


E 


H 


X 


SU(2) x U(l) 


2i/6 I2/3 


1-1/3 


2-1/2 


l-i 


2l/2 


lo 


A 4 


singlets singlets 


singlets 


3 


3 


(MM") 


1 



Table 2: The A4 charge assignment of the SM fields and the dark matter X. 



the four-Fermi decay interaction is explicitly written as 

•^decay 



(LE)LX + ^ (L£)'LX + h.c. 



(2.2) 



A. 



aT ( v r^R ± ^Tb) e L X - (t l /j, r ± /i L r fl ) z/ e X 



+ (v 6 tr ± u T e R ) fM L X - (e L r R ± r L e R ) v^X 

+ (wjle R ± V^fi R ) TEX - (jZEe R ± eZ^s) T\X + h.c. (2.3) 

There are two types of operators, which we have denoted with the coefficients A±, 
corresponding to the fact that there are two ways to construct the A 4 triplet represen- 
tation from two 3's. It should be noted that, due to the non-Abelian A 4 symmetry, 
the decay vertices have specific structures of chirality and generations. 

We have introduced three higgs doublets H\ y in to have the masses of charged 
leptons (the details of lepton masses and mixing will be discussed in later section). 
It was shown [T2] that the introduction of multi higgs doublets in this manner does 
not lead to dangerous flavor- changing processes. 



3 Cosmic-ray anomaly 

In this section, we show by calculating the positron (electron) flux that the scenario 
given above, which has a special generation structure of DM decay vertices, is possi- 
ble to excellently describe the cosmic-ray anomalies reported by the PAMELA and 
Fermi-LAT experiments. 



3.1 Positron production from DM decay 

First, we consider the branching fraction of the DM decay through the A 4 - invariant 
operator LELX. Due to the typical generation structure given in ( 12. 3p . the dark 
matter X decays into several tri-leptons final state with the equal rate, where each 
final states include all three flavors: 

Br(X -> e ± ^v T ) = Br(X -> r+e^) = Br(X -> ^T^v e ) = -. (3.1) 
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Here we have omitted the masses of charged leptons in the final states. The branching 
fractions indicate that the spectrum of positrons (electrons) in cosmic rays is uniquely 
determined in the present framework with A4 symmetry, which allows us to predict 
the spectrum of cosmic-ray anomalies. The total decay width of DM turns out to be 

where mx is the DM mass. 

Given the decay width and the branching fractions, the positron (electron) pro- 
duction rate (per unit volume and unit time) at the position x of the halo associated 
with our galaxy is evaluated as 



Q(E,x) = n x (x)Tj2MX -> /) 

/ 



dE 



(3.3) 
J/ 



where [dN e ±/dE]f is the energetic distribution of positrons (electrons) from the decay 
of single DM with the final state '/'. We use the PYTHIA code [13J to evaluate the 
distribution [dN e ±/dE]f. The DM number density n x (%) is obtained by the profile 
p(x), the DM mass distribution in our galaxy, through the relation p(x) = mxnx{x). 
In this work we adopt the Navarro-Frank- White profile 



/ -*\ r Q (r & + r c ) 2 

PNFW{X) = p Q 1— rg-, (3.4) 

r(r + r c ) A 

where p & ~ 0.30 GeV/cm 3 is the local halo density around the solar system, r is the 
distance from the galactic center whose special values r ~ 8.5 kpc and r c ~ 20 kpc 
are the distance to the solar system and the core radius of the profile, respectively. 

In the present model, the dark matter decays into not only and /j, ± which 
result in pure leptonic decays, but also r 1 * 1 leading to hadronic decays, and anti- 
protons may also be produced in the halo of our galaxy. It is however obvious that 
the dominant decay channels are leptonic and the branching fractions of hadronic 
decay are made tiny by the electroweak coupling and the phase space factor. The 
suppression of hadronic decays is consistent with the p data obtained in the PAMELA 
experiment [5]. On the other hand, the injections of high-energy positrons (electrons) 
in the halo give rise to gamma rays through the bremsstrahlung and inverse Compton 
scattering processes. Comprehensive analyses of cosmic-ray fluxes [5] show that the 
gamma-ray flux from leptonically decaying DM is also consistent with the Fermi- 
LAT data [13] . As a result, we concentrate on the calculation of positron (electron) 
flux in what follows. 
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3.2 Diffusion model 

Next, we consider the propagation of positrons (electrons) produced by the DM 
decay in our galaxy. The charged particles e suffer from the influence of tangled 
magnetic fields in the galaxy before arriving at the solar system. The physics of the 
propagation can be described by the diffusion equation [T61 E] > 

K e± (E)V 2 f e ±(E,x) + -^[b(E)f e± (E,x)} + Q(E,x) = 0. (3.5) 

The number density of per unit energy f e ± , satisfies the condition f e ± = at the 
boundary of the diffusion zone. The diffusion zone is approximated to be a cylinder 
with the half-height of 4 kpc and the radius of 20 kpc. The diffusion coefficient 
K e ±(E) and the energy-loss rate b(E) are set to be 

K e ±(E) = 1.12 x 1(T 2 [kpc 2 /Myr] x £° e 7 v °, (3.6) 
b(E) = 1.00 x 10~ 16 [GeV/sec] x #| eV , (3.7) 

where E Ge y = E /(l GeV). To fix these parameters, we have used the MED set for 
the propagation model of [18J, which gives the best fit value in the boron-to- 
carbon ratio (B/C) analysis as well as in the diffused gamma-ray background. Once 
f e ± is determined by solving the above equation, the e 1 * 1 fluxes are given by 

[<M£)]dm = ^f e ±(E,x Q ), (3.8) 

where x & is the location of the solar system, and c is the speed of light. For the 
total fluxes of e ± 1 we have to estimate the background fluxes produced by collisions 
between primary protons and interstellar medium in our galaxy. In the analysis, the 
following fluxes for cosmic-ray electrons and positrons [J7] are adopted: 

r$ 1 . = O- 16 ^ 1 (39) 

1 -f ll-C/ GcV -+- o.z_C/ GeV 



[$ i = °- 70 ^c 7 v (3 10) 

^e-j sec 1 + 110^ + 600^ + 580^' {6 Uj 

A.5E - 7 

[$e+]sec = 1 + 650^ + 1500^' (3 ' n) 

in unit of (GeV cm 2 sec str) -1 . With these backgrounds, the total fluxes and the 
positron fraction R e +, which is measured by the PAMELA experiment, are found to 



6 




Figure 1: The positron fraction and the total e + + e~ flux predicted in the leptonically- 
decaying DM scenario with A4 symmetry. The DM mass is fixed to 1, 1.5, and 2 TeV. As 
for the DM decay width used in the fit, see the text. 



be 

[$e+]total = [$e+]DM + [$e+]sec, (3.12) 

[$e-]total = [$e+]DM + a[$e-]prim+ [$e-]sec, (3.13) 

R e + = [$ e +]total/([$ e+] total + [$e-]total). (3.14) 

Note that the primary flux for electrons measured by Fermi-LAT should be multiplied 
by the normalization factor a = 0.7 so that our evaluation is consistent with the 
experimental data in the low-energy range |19j . 



3.3 Results for PAMELA and Fermi-LAT 

The positron fraction and the total flux [3> e -]totai + [3> e +] total are depicted in Figure [T] 
for the scenario of the leptonically decaying DM with A4 symmetry. For the DM 
mass nix = 1, 1-5, and 2 TeV, the results are shown with the experimental data of 
PAMELA and Fermi-LAT. The total decay width T is fixed for each value of DM 
mass so that the best fit value explains the experimental data. With a simple x 2 
analysis, we obtain T _1 = 1.7 x 10 26 , 1.2 x 10 26 , and 9.5 x 10 25 sec for m x = 1, 1.5, 
and 2 TeV, respectively. It can be seen from the figure that the PAMELA anomaly 
is well explained in the decaying DM scenario with A 4 symmetry. Furthermore, the 
latest Fermi-LAT data is perfectly fitted in this scenario if the DM mass is around 
2 TeV. 
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4 Lepton masses and mixing 



So far, the dark matter property, especially the leptonic decay, has been analyzed 
for the gauge-singlet fermion X. In this section, we discuss the lepton masses and 
mixing in the same setup as Table [2] and also in two types of its extensions. 

For the matter content and the A4 assignment given in Table EJ the charged- 
lepton and neutrino masses come from the symmetry-invariant operators 



C 



£ (y e )iH t LE + h.c.+ {Vvh LcH t*HfL. 

i,j=l,V,l" 



(4-1) 



8=1,1' 1 



The subscripts i mean the singlet representations of A 4 symmetry, i = 1, V, 1". The 
higgs fields are assumed to develop vacuum expectation values (Hi) = (0, Uj/\/2) t - 
The lepton mass matrices turn out to take the forms 



M P 



m e 



\ 



( 



mi 



m, 



m 2 



m. 



\ 



m 3 J 



(4.2) 



m fj. = f{vu uvy, lu 2 Vi»), m 2 = g(vi,Lo 2 vii,Lovin), 

m r = f(vi,U! 2 Vi',U!Vi"), ins = g(vi,cuvi',uj 2 Vi>>), (4.3) 

where to = e 2m / 3 , and the functions / and g are given by 

1 . 1 . 

f(v u v v ,v v ,) = -T=22{y e )iVi, g(vuv v ,v v ,) = 2z2(y»h v i v r ( 4 - 4 ) 

With suitable values of the coupling constants, the experimentally-observed masses 
(differences) are able to be reproduced^ The generation mixing is, however, ab- 
sent unless some ingredient is added. In the following, we will present two possible 
examples to remedy this problem without causing a rapid decay of the dark matter. 

The first example is to introduce extra higgs doublets which induce Majorana 
neutrino mass, i.e., additional dimension 5 operator like (14. li . The extra higgses H' 
belong to the triplet representation of A 4 symmetry in order for non-trivial flavor 
mixing to be generated. Further, H should be charged under some symmetry not 

When Vi are the electroweak scale, the neutrino mass M v - lCT (1 - 2) eV seems to imply that 
the effective scale of LLHH operator, y~ x ~ A', is somewhat below the unification scale, namely, 
the lepton number symmetry is valid above A' in low-energy effective theory. 
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to have the interactions (the effective operators listed in Table [T]) which cause the 
DM decay and disturb the previous result. To satisfy this requirement, we consider 
a simple example with Z2 parity under which only H' is negative. As a result, the 
decay operators involving H' with dimensions less than 7 are not permitted, except 
for dimension 6 operators H'* H'LH C X and H'D^H'X^^E. It is found that they 
cannot be forbidden by any Abelian (discrete) symmetry while other necessary terms 
remain intact. Therefore, if one assumes that the DM decay from these operators is 
sub-dominant, the expectation values of H' should be suppressed. 

The remaining is the decay operator of dark matter LELX and the additional 
source of neutrino masses y' u L c LH'H'. The charged-lepton masses are unchanged 
and the neutrino mass matrix turns out to be 



M' 



( m[ y' u v' 2 v' 3 y' u v[v' 3 \ 
y' v v' 2 v' z m' 2 y' v v[v' 2 



(4.5) 



\y'A v 3 m 3 J 

where v[ are the expectation values of H[. The diagonal elements vn! i are shifted by 
0(y' v vf) from rrii due to the new interaction, and their exact forms are determined 
by the A4 invariance. The additional 3 degrees of freedom (the off-diagonal matrix 
elements) can fit the experimental values of neutrino mixing. 

Another way to have non-vanishing generation mixing is to consider a different 
type of neutrino mass operator than (14. Xp with use of the SU(2)-triplet scalar A. 
The simplest tree-level term for neutrino mass is constructed with the scalar A: 

C A = yJs^L. (4.6) 

Similar to the first example, A should belong to the triplet representation of A4 
symmetry for non-trivial generation mixing of neutrinos. The electroweak gauge 
invariance implies that the above term only induces off-diagonal elements in the 
neutrino mass matrix. Assuming nonzero expectation values va% for the neutral 
components of A, (i = 1, 2, 3), we obtain the neutrino mass matrix 



/ mi y A VA3 y^VA2^ 



M A 



(4.7) 



y&VA3 m 2 l/A^Al 
\yAVA2 VaVai ^3 ) 

The phenomenological analysis based on this type of Majorana mass matrix has been 
performed in Ref . [20] , where the solar and atmospheric neutrino anomalies and the 
neutrino-less double beta decay have been studied. 
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It is noticed that the triplet scalar A gives rise to new decay interactions of dark 
matter. For operators with dimensions more than 5, their contributions to the decay 
amplitude are suppressed when A is heavier than the dark matter and va is much 
smaller than Vi to satisfy the electroweak precision (the p parameter constraint). 
The gauge and flavor invariance then leave a single dimension 5 operator 

\ A H£*LX. (4.8) 

It is easily found that this operator cannot be forbidden by imposing any symmetry, 
if one allows the necessary operators for the lepton masses and the DM decay through 
LELX. To avoid a rapid DM decay via the operator ( 14. 8|) . v& should be smaller than 
(TeV) 2 /A ~ eV. Then the coupling y& in (14.61) is 0(1) for non-negligible neutrino 
mixing. Integrating out the heavy scalar A with its mass m^, we have an effective 
operator 

Va A a 



Am A 



HL C LLX. (4.9) 



Since the coupling y& is 0(1), this dimension 7 operator might give a sizable effect 
on the X decay. In other words, if one requires that the dominant decay vertex is the 
four-Fermi operator LELX, the triplet scalar should be heavier than the intermedi- 



ate scale: wa > y |Aa|Au . Such an SU(2)-triplet scalar with an intermediate mass 
and a tiny expectation value might be incorporated in SO(10) unified theory with 
the intermediate Pati-Salam group, where the potential analysis is slightly shifted 
by the electroweak scale. We finally mention that a tiny value of Aa (< TeV/A) 
might also be a solution with low-mass A. That however means the effective theory 
description is invalid and the model should be improved. 



5 Conclusion 

We have considered the decay of gauge-singlet dark matter for the cosmic-ray anoma- 
lies reported by the PAMELA and Fermi-LAT experiments. The decaying dark mat- 
ter recently attracts much attention because, if it is a TeV-scale fermionic particle, 
the suggested order of meta-stability is just derived from a four- Fermi interaction 
suppressed by the GUT scale. It is also noted that the cosmic-ray anomalies are 
explained by the DM decay, while the relic abundance may be determined by DM 
annihilation process, e.g. mediated by a light singlet scalar. 
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The scenario is however spoiled due to the existence of other operators which 
force a rapid DM decay and/or induce non-leptonic DM decay. In this letter, we 
have pointed out that such harmful decay vertices are prohibited by implementing 
a single non-Abelian flavor symmetry such as A 4 . Any Abelian symmetry cannot 
play the same role. We have also shown that the A 4 invariance leads to the flavor- 
universal decay channels of DM, with which the cosmic-ray anomalies are captured 
very well with the DM mass around 2 TeV. Further we have discussed the relevance 
of discrete flavor symmetry on neutrino phenomenology and offered two independent 
mechanisms to generate lepton masses and mixing without disturbing the successful 
decaying DM scenario. It would be therefore interesting to construct a high-energy 
completion, i.e. a concrete GUT model involving both the dark matter candidate 
and mechanism to generate neutrino masses with a non-Abelian discrete symmetry. 
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